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Upsilon Decays into Scalar Dark Matter 
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We examine decays of a spin-1 bottomonium into a pair of light scalar Dark Matter (DM) particles, assuming 
that Dark Matter is produced due to exchange of heavy degrees of freedom. We perform a model-independent 
analysis and derive formulae for the branching ratios of these decays. We confront our calculation results 
with the experimental data. We show that Dark Matter production in T decays may lead to constraints on 
parameters of the models containing a light spin-0 DM particle. 



1. Introduction 

We consider the possibility of using of T meson 
decays with missing energy, to test the models with 
a light spin-0 DM particle. Dark Matter search in 
heavy meson (and in particular T meson) decays 
may be complimentary to such experiments as DAM A 
[3,0,3, CDMS [| and XENON 0,0, which rely on 
the measurement of kinematic recoil of nuclei in DM 
interactions and lose (for cold DM particles) sensitiv- 
ity with decreasing mass of the WIMP, as the recoil 
energy becomes small. 

So far, T meson decays into Dark Matter have been 
considered within the models, where DM particles in- 
teraction with an ordinary matter is mediated by some 
light degree of freedom @, [8, 9]. Apart from desire of 
having DM annihilation enhancement (due to a light 
intermediate resonance) and thus having no tension 
with the DM relic abundance condition [l(|-[l3j], it is 
also known that T meson SM decay is predominantly 
due to strong interactions. Thus, the WIMP produc- 
tion branching ratio, in general, is greatly suppressed 
compared to relevant weak B decays, and in particu- 
lar to B — > K + invisible transition [Til . [l5| . In light 
of this, it might seem natural to concentrate only on 
the models within which Dark Matter production in T 
decays is enhanced due to exchange of a light particle 
propagator. 

Yet, our aim is to study T(IS') decay into a pair of 
spin-0 DM particles, T(15) -> and T(35) de- 

cay into a pair of spin-0 DM particles and a photon, 
T(3S) — > $$*7, within the models where light Dark 
Matter interaction with an ordinary matter is due to 
exchange of heavy particles (with masses exceeding 
the bottomonium mass). These models may be free of 
tension related to satisfying the DM relic abundance 
constraint as well [HI E3, EE E3, EH- Also, new ex- 
perimental data on T decays into invisible states have 
been reported by the BaBaR collaboration (l9l . [20j . 
According to these data, 

S(T(15) -> invisible) < 3 x 10~ 4 (1.1) 

and 

B(T(3S) -> 7 + invisible) < (0.7 - 31) x 10~ 6 (1.2) 



where the interval in the r.h.s. of eq. (| 1 . 2[) is related 
to the choice of the final state missing mass. These 
bounds are significantly stronger than those on invisi- 
ble T(15) decays (with or without a photon emission), 
reported previously by Belle and CLEO [U [22j and 
quoted by Particle Data Group [23j. We show that 
BaBaR experimental data on T meson invisible de- 
cays may constrain the parameter space of light scalar 
Dark Matter models, even if there is no Dark Mat- 
ter production enhancement due to light intermediate 
states. 

We also illustrate that the study of Dark Matter 
production in T decays allows us to test regions of 
parameter space of light spin-0 DM models that are 
inaccessible for B meson decays with missing energy. 
It is also worth mentioning that T decays are sensi- 
tive to a wider range of WIMP mass than B decays. 
Thus, the study of WIMP production in T decays is 
complementary to that for B meson decays. 

At the energy scales, associated with T decays, 
heavy intermediate degrees of freedom may be inte- 
grated out, thus leading to a low-energy effective the- 
ory of four-particle interactions. Our strategy would 
be deriving first model-independent formulae for the 
T(15) -> and T(35) -> $$*7 branching ra- 

tios within the low-energy effective theory. Then, we 
confront our predictions with the experimental data, 
deriving model-independent bounds in terms of the 
Wilson operator expansion coefficients, as the param- 
eters that carry the information on an underlying New 
Physics model. Finally, within a given model, using 
the matching conditions for the Wilson coefficients, 
we translate these bounds into those on the relevant 
parameters of the considered model. 

The talk is based on the results presented in [24[ 
and [U. 



2. Model-Independent Analysis 

We treat T states - neglecting the sea quark and 
gluon distributions - as bound states of bb valence 
quark- antiquark pair that annihilates - with or with- 
out emission of a photon - into a pair of Dark Matter 
particles. To this approximation, the relevant low- 
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Figure 1: Upper bound on [C3 ] a) as a function of m$, for Ah = 100 GeV, b) as a function of Ah, for m$ = 1 GeV. 



energy effective Hamiltonian may be written as 



(2.1) 



where Ah is the heavy mass and 
Oi = m b (6 b) ($*$) , 2 = im b (b l5 b) ($**) , 
3 = (frfb) (**ifl p $) , 4 = (67^756) (V^$) (2.2) 

with 9 = 1/2(9 — 9). It is worth noting that with the 
notations used in (|2.ip and (|2.2p . all the operators Oi, 
i=l,..4, are Hermitean, thus all the Wilson coefficients 
Ci must be real. If DM consists of particles that are 
their own antiparticles, then only first two operators 
in eq. (|2.2p would contribute. 

For DM field being a complex scalar state, T(15) — > 
$$* branching ratio is given by the following expres- 



B(T(1S) 



$$*) = 



r(T(lS) -► $$*) 



T(1S) 



,f 2 



T(1S) 



A|j 487rr T 



(1S) 



M- 



T(1S) 



— 4m 



.3/2 



(2.3) 



where 771$ is the DM particle mass and r T (is)j 
-Wr(is), /t(is) are the T(15) total width, mass and 
decay constant respectively. Only operator O3 con- 
tributes to the decay rate. 

For $ being a self-conjugate spin-0 state, $ = $*, 
contribution of O3 vanishes as well, as noted above. 
Thus, one has 



S(T(15) -> $$) = 



(2.4) 



This result is related to the fact that the final DM par- 
ticle pair state must be a P-wave, which is impossible 
due to the Bose-Einstein symmetry of identical spin- 
particles. In what follows, r(T(15) -> $$) must 
also vanish in higher orders in 1/ mj, operator product 
expansion. 



Thus, provided that DM pair production is the 
dominant invisible channel (the neutrino background 
may be neglected [M III III ) , the signal for T(15) ->■ 
invisible decay would imply that the light spin-0 DM 
field has a complex nature. No evidence for the 
T(15) — > invisible mode may lead to some constraints 
on the parameters of the models with light complex 
scalar Dark Matter. 

Yet, in order to derive such constraints, the experi- 
mental limit on the T(1S) — > invisible mode must be 
strong enough. Indeed, using the numerical values of 



^t( is); 1 Mr (is) 
eq. (H3D as 



and f 



T(1S) 



[26l |. one may rewrite 



B(T(IS) -> $$" 
ZlOOGeV 
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(2.5) 



In what follows, the relevant experiments must be sen- 
sitive (at least) to B(T(1S) — > invisible) ~ 10 -4 . 

This sensitivity has been reached by the BaBaR 
experiment [T^ |. as it follows from the bound on 
B(T(1S) — > invisible), given by eq. (|1.1[) . Substi- 
tuting (|1.1[) into (|2.5p . one derives the following con- 
straint on IC3I as a function of m$ and Ajy: 



IC3I < 0.75 



A 



(ioo<rv< 1 



4m| 
M 2 

iw T(lS) , 



-3/4 



(2.6) 



The behavior of the upper bound on \Cs\ with the 
DM particle mass and the heavy mass is presented in 
Fig.'s[T^) andQJ)) respectively. For m$ < 3 GeV and 
Ah ~ 100 GeT^, this bound may be translated into 
constraints on the relevant couplings of models with 
a light complex spin-0 DM field. For to$ > 3 GeV or 
for Ah > 100 GeV, further improvement of the exper- 
imental sensitivity to T(1S) — > invisible transition is 
necessary. 

Note that even for m$ < 3 GeV and Ah — 
100 GeV, bound still allows C 3 to be of the 
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order of unity. It seems to be very unlikely to sat- 
urate such a (rather weak) bound, if within a given 
New Physics model, T(15 l ) — * transition is loop- 
induced. Thus, bound (|2.6|) on | C3 1 seems to be useful 
only within the models, where T(IS') — > $$* decay 
occurs at tree level. 

Further discussion of the scenarios with complex 
scalar Dark Matter goes beyond the scope of this talk. 
Yet, in ref. [24[ we use bound (|2.6[) to derive con- 
straints on the parameter space of the models with 
mirror fermions, in particular on that of the MSSM 
with gauge mediated SUSY breaking and the DM par- 
ticle in the hidden sector, while mirror fermions being 
connectors between the hidden and the MSSM sectors. 
We derive for the first time bounds on the couplings 
of b-quark interactions with its mirror counterpart, as 
functions of the WIMP mass and the mirror fermion 
mass. 

For the scenarios with the DM particle being its own 
antiparticle, the decay T(35 l ) — > $$7 is relevant. As 
the neutrino background is negligible [24J, T(35) — * 
$$7 may be the dominant channel in the T(35) — * 
7 + invisible mode. 

As it was mentioned above, the strongest experi- 
mental constraint on B(T — ► 7 + invisible) is that 
derived by the BaBaR collaboration [2(| and given by 
cq. (|1.2p . It may also be written as 



A H = 100 GeV 



B(T(3S) -> 7 + invisible) < 3 x 10" 



(2.7) 



for the final state missing mass being less than 7 GeV, 
or approximately ^/s < Mx(3S)/v2- We have to 
note however that this bound has been derived as- 
suming that T(3S) — > 7 + invisible transition is me- 
diated by an intermediate resonant Higgs state A 
that exists in the MSSM extensions with an addi- 
tional Higgs singlet 0, EH, H3] ■ For the models with 
non-resonant DM production considered here, bound 
(|2.7p may be used only to get a preliminary estimate 
of possible constraints on parameters of light spin- 
self-conjugate DM models 1 . More rigourously, the 
experimental analysis performed in ref. [2(| should 
be extended to the cases, when the emitted pho- 
ton energy is non-monochromatic and is in the range 
< uj < Mx(3S)/2 — 2m|/Af T ( 3 5). To our knowledge, 
this work is in progress now 2 . 

Within self-conjugate scalar DM scenarios, the par- 
tially integrated branching ratio for T(35) — > $$7 



1 As for the earlier bounds on BIT — > 7 + invisible) [23ll . 
the only existing constraint for the case of two-particle invis- 
ible states and the photon having non-monochromatic energy, 
B(T(15) -> 7+11) < 10~ 3 [53, is too weak to have any use. 
Due to the factor a/(47r) w 5.8 x 10~ 4 , B(T -> $$7) is always 
below this limit. 

2 Author is grateful to Yu. Kolomensky for this information. 
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Figure 2: Upper bound on \J C\ + C| as a function of m$, 
for A H = 100 GeV. 



decay is given by 



(C? + Cf) a I\ 



(3S) 



A H 47T 547rr T (35)M T (35) 



^T(3S) 




2m$ 

For s maa; = M^ 3S -j/2, using the numerical values of 
the T(35) mass, total width and decay constant [23|, 
[24[, one may rewrite eq. (|2.8p in the following form: 

B(T(3S) -» $$7) 
= 2.6 x 



(2- 



I S < M T(3S)/ 2 



ID" 7 (C 2 + C 2 ) (^g^) /(«.) (2-9) 



where 



and 



/(a*) = ( 1 + -£* ) Vl - 8x$ 



" ygft(l-g»)ln ^ 2 ^^_ j (2.10) 

At first glance, it may seem that T(35) — > $$7 
branching ratio is far out of reach of the BaBaR ex- 
perimental sensitivity for a reasonable choice of C\ 
and C2. Notice, however, that within certain mod- 
els with light spin-0 Dark Matter, the Wilson coeffi- 
cients C\ and/or Ci may be enormously large, as they 
contain some enhancement factors, such as the ratio 
Ah I TOj 3> 1 - due to the mass term in the numerator 
of a heavy fermion propagator, or the Higgs vev's ra- 
tio tan/3 (with the latter being, say, ~ mt/mb 3> 1) - 
due to DM particle pair production via exchange of a 
heavy non-SM Higgs degree of freedom. 

These enhancement factors can make C\ and/or C2 
to be > 10 and hence B(T(3S) -> $$7) to be ~ 
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Figure 3: Upper bound on |A h o s | as a function of m$, for 
m h o = 55 GeV and tan/3 = 20 (line 1), tan/3 = 40 (line 
2). 



10 5 — 10 , i.e. significantly exceeding bound (|2 . T[) 
on B(T — » 7 4- invisible). Using bound (|2.7[) yields 

The behavior of the bound on ^/C\ + Cf with the DM 
particle mass is presented in Fig.[2]for Ah = 100 GeV. 

In the next section we consider the simplest pos- 
sible model where the Wilson coefficients are made 
enormously large due to an en hancemen t factor. We 
transform the constraint on \J C\ + C\ into that on 
the relevant parameters of the model. 



3. Dark Matter Model with two Higgs 
Doublets (2HDM) 

In this section we consider the two-Higgs doublet 
model (2HDM) with a gauge singlet real scalar DM 
particle. The DM interaction part of Lagrangian, rel- 
evant for our analysis, may be written as [14j 

_ C = 4- Ar^ltfxl 2 4- A 2 $ 2 |H 2 | 2 + 

4-A 3 $ 2 (HiH 2 + h.c) (3.1) 

Following refs. 0, H8[, we consider type-II version 
of 2HDM, where H\ generates masses of down-type 
quarks and charged leptons, whereas H 2 generates 
masses of up-type quarks. We assume for the Higgs 
vev's ratio to be large, i.e. Vv/vx = tan/3 3> 1. 

Three coupling constants entering Lagrangian (|3.ip 
are unknown parameters. Constraints on the cou- 
plings Ai and A 2 may be derived from the study of 
B — > K 4- invisible transition. Those, in general, are 
strong enough: combined with the ones coming from 
the DM relic abundance condition, they rule out a 
wide range of the WIMP mass for the scenarios with 
Ai and/or A 2 dominant (l4| . 



Yet, due to cancellation effects in the relevant di- 
agrams, WIMP pair production rate in B meson de- 
cays is insensitive to the value of A3 [3]. The sce- 
narios with A3 dominant, or at least non-negligible, 
are thus far unconstrained. Study of DM production 
in T decays enables one to constrain these scenarios, 
inaccessible by B meson decays with missing energy. 

The model contains two CP-even, one CP-odd and 
two complex charged physical Higgs degrees of free- 
dom. To the leading order in the perturbation the- 
ory, Y(35) — * $$7 transition occurs at tree level by 
exchange of a single Higgs boson. To this approxima- 
tion, only the CP-even Higgs states, h° and H°, are 
relevant for our analysis. 

One of the CP-even Higgs bosons (not-necessarily 
the lightest one) is Standard Model-like: its phe- 
nomenology is similar to that of the SM Higgs boson 
and the experimental bound on its mass is close to the 
SM limit 3 (see [23[ and references therein) . Contribu- 
tion of the diagrams with the SM-like Higgs boson 
exchange to the T(3S) — > $$7 decay amplitude has 
no any enhancement factor - thus, we may further dis- 
regard the Dark Matter interaction with the SM-like 
CP-even Higgs. 

The other CP-even Higgs boson is "New-Physics 
(NP) like": its phenomenology differs drastically from 
that of the SM Higgs boson [31]. If Dark Matter 
is produced due to exchange of this Higgs particle, 
T(35) — * $$7 decay amplitude is enhanced by tan/3 
factor. As discussed above, in this case the decay 
branching ratio is within the reach of the BaBaR ex- 
perimental sensitivity. 

If the NP-like Higgs is the lightest one, its mass 
may be much below the Standard Model experimen- 
tal limit: according the existing experimental data 
[H, HU, m h o > 55 GeV or m h o < 1 GeV in the gen- 
eral type II 2HDM. As it was mentioned above, the 
light Higgs scenario is beyond the scope of the present 
paper, thus we assume here that m^o > 55 GeV. No- 
tice, however, that this bound is derived, provided 
that no invisible Higgs decay mode exists. On the 
other hand, if the NP-like Higgs invisible decay mode 
is dominant, it may escape detection. No bound on 
m/jo, to our best knowledge, exists in that case. 

Within the considered model with light scalar Dark 
Matter, analysis of the T — ► $$7 mode may restrict 
the scenarios with an invisibly decaying lightest Higgs 
boson by putting severe constraints on the in- 
teraction coupling, AfjO$. For tan/3 ^ 1 it is given 
by 

Ko<s> ~ A 3 4- (Ai 4- A 2 )cos/3 (3.2) 



3 Also, if the SM Higgs decays predominantly invisibly, the 
SM lower experimental bound is distorted by a few GeV only 

mm. 
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a) m H D = 160 GeV 
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Figure 4: Upper bound on | A3 1, a) as a function of m$ for m H o — 160 GeV and tan — 30 (line 1), tan /3 = 55 (line 2), b) 
as a function of m H o for tan/3 = 55 and m<j = 100 MeV (solid line), m$ = 1.5 GeV (dashed-dotted line), m<j = 2 GeV 
(dashed line). 



The last term in the r.h.s. of (|3.2j) . although being 
suppressed by a factor of cos/3 « 1/ tan/3, must be 
retained because of possible hierarchy in the values of 
A3 and Ai or A2. Scenarios with such a hierarchy may 
be of importance: using the matching conditions for 
the Wilson coefficients, 

C\ = --^jptan/3, C 2 = 0, A H = m h o, (3.3) 



and bound (|2.1ip on \/Cf + Cf , one finds that A^o$ is 
constrained to be 0(1/ tan/3), if h° mass approaches 
to its lower limit, m^o = 55 GeV. We illustrate this 
in Fig. rj for tan/3 = 20 and tan/3 = 40. 

For A/jO$ being so strongly constrained, it seems to 
be very unlikely that h° would escape detection and 
its mass be below 55 GeV. More rigorously, however, 
detailed re-analysis of the Higgs production and de- 
cay rates, including that of h° — > $$, should be per- 
formed. 

One can show [24| that when choosing tan/3 suffi- 
ciently large, bound on A^o$ may still be rigorous, if 
the ft, mass is heavier than 55 GeV. Thus, within the 
type II 2HDM with a light spin-0 Dark Matter, study 
of T — > $$7 decay channel may lead to severe con- 
straints on the lightest CP-even Higgs invisible decay 
coupling, if that Higgs is New-Physics like. 

The remarkable feature of the model is that con- 
straints on the parameter space are derived even if 
the NP-like Higgs is the heaviest CP-even one. We 
illustrate this choosing m H o > 160GeV - the exist- 
ing theoretical upper bound on the SM-like Higgs bo- 
son [34| and exclusion of the SM Higgs mass interval 
(160- 170) GeV with 95% C. L. by the CDF and DO 
data |35j allow us to infer that above 160 GeV, the 
CP-even Higgs boson is presumably the heaviest one 
and NP-like. The matching conditions for the Wilson 
coefficients are now 

Ci = -Vtan/3^ C2 = 0; AH = mm (34) 



Thus, bound (|2~TTj) on y/C\ + Cf may be transformed 
into that on | A3 1 as a function of the WIMP mass, 
tan/3 and the heaviest CP-even Higgs mass. 

As one can see from Fig. [4^,), for m#o = 160 GeV 
and tan/3 = 30, A3 is constrained to be of order of the 
SM weak coupling or smaller (| A3 1 < 0.65), if m$ < 
1 GeV . Also, for the same choice of the Higgs mass 
and tan/3, | A3 1 is to be less than one, if the WIMP 
mass is less than 2 GeV. Bound on | A3 1 is significantly 
more rigorous for higher values of tan /?. For instance, 
if choosing tan/3 = mt(mt) / mb(m t ) » 55, one gets 
|A 3 | < 0.35 and |A 3 | < 0.5 for to* < 1 GeV and 
m$ ~ 2 GeV respectively. 

The restrictions on | A3 1 are essential also for higher 
values of the heaviest CP-even Higgs mass: for tan (3 = 
55, they are still of the interest up to m H o ~ 280 GeV, 
as one can see from Fig. |4}d). 

Thus, within the type II 2HDM with a scalar Dark 
Matter, for large tan (3 scenario, T meson decay into a 
Dark Matter particles pair and a photon, T — > $$7, 
may be used to derive essential constraints on the 
parameters of the model, which otherwise cannot be 
tested by B meson decays with invisible outcoming 
particles. 



4. Conclusions and Summary 

Thus, spin-0 Dark Matter production in T meson 
decays has been investigated. We restricted ourselves 
by consideration of the models where the decays oc- 
cur due to exchange of heavy non-resonant degrees of 
freedom. 

We performed our calculations within low-energy 
effective theory, integrating out heavy degrees of free- 
dom. This way we derived model-independent for- 
mulae for the considered branching ratios. We used 
these formulae to confront our theoretical predictions 
with existing experimental data on invisible T decays, 
both in a model-independent way and within particu- 
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lar models. It has been shown that within the consid- 
ered class of models, DM production rate in T decays 
is within the reach of the present experimental sensi- 
tivity. Thus, T meson decays into Dark Matter, with 
or without a photon emission, may be used to con- 
strain the models with a GeV or lighter spin-0 DM. 

Experimental constraints on the T(3S I ) — > 7 + 
invisible mode are derived assuming that Dark Mat- 
ter is produced by exchange of a light resonant scalar 
state. Within the scenarios with non-resonant DM 
production considered here, these constraints may be 
used only to make preliminary estimates of possible 
bounds on the parameters of the models. Yet, those 
estimates show that these bounds may be rigorous 
enough; besides, they are derived within the least re- 
strained presently light scalar DM scenarios. Our goal 
is thus to encourage the experimental groups to ana- 
lyze the experimental data on T — > 7 + invisible also 
for the case of non-monochromatic photon emission 
and spin-0 invisible states. 

So, from our analysis one may conclude that Dark 
Matter production in T meson decays may serve as an 
interesting alternative to commonly used DM search 
methods, capable of providing a valuable information 
on DM particles, if those turn to have a mass of the 
order of a few GeV or smaller. 
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